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ABSTRACT 

We present a detailed analysis of the "second-parameter pair" of globular clusters M3 (NGC 5272) and 
Palomar 3. Our main results can be summarized as follows: i) The horizontal-branch (HB) morphology 
of M3 is significantly bluer in its inner regions (observed with the Hubble Space Telescope) than in the 
cluster outskirts (observed from the ground), i.e., M3 has an internal second parameter. Most plausibly 
the mass loss on the red giant branch (RGB) has been more efficient in the inner than in the outer 
regions of the cluster, ii) The dispersion in mass of the Pal 3 HB is found to be very small — consistent 
with zero — and we argue that this is unlikely to be due to a statistical fluctuation. It is this small mass 
dispersion that leads to the most apparent difference in the HB morphologies of M3 and Pal 3. iii) The 
relative HB types of M3 and Pal 3, as measured by mean colors or parameters involving the number of 
blue, variable, and red HB stars, can easily be accounted for by a fairly small difference in age between 
these clusters, of order 0.5-1 Gyr — which is in good agreement with the relative age measurement, based 
on the clusters' turnoffs, by VandenBerg (2000). 

Subject headings: Hertzsprung-Russell (HR) diagram and C-M diagrams — stars: horizontal-branch — 
stars: mass loss — stars: Population II — globular clusters: individual: Palomar 3, 
M3 (NGC 5272) 



1. INTRODUCTION 

Many recent studies of the globular clusters (GCs) in 
the outer halo of our Galaxy have focused on the age dif- 
ference between these distant clusters and those located 
in the inner halo. Since Mironov & Samus (1974) and 
Searle & Zinn (1978), it has been widely recognized that 
the possible existence of such age differences would have 
important implications for our understanding of the early 
stages of the formation of the Galaxy (e.g., van den Bergh 
1993; Zinn 1993; Lee, Demarque, & Zinn 1994; Majewski 
1994). While in these early studies the morphology of the 
horizontal branch (HB) was assumed to be a safe relative 
age indicator — a critical working hypothesis in those days, 
given the tremendous challenge that accurate and precise 
photometry reaching below distant clusters' turnoffs (the 
primary stellar evolution "clock") represented — more re- 
cently Hubble Space Telescope (HST) photometry of the 
most distant GCs in the Galaxy have allowed a direct 
and precise estimate of the position of their main-sequence 
turnoffs, thus enabling the direct determination of their 
ages as compared to closer GCs with good ground-based 
photometry (Harris et al. 1997; Stetson et al. 1999; Van- 
denBerg 2000). 

As a result, it was found that the most massive of the ex- 
treme outer-halo GCs, NGC 2419, has the same age as the 
well-studied metal-poor, inner-halo GC M92 = NGC 6341 
(Harris et al. 1997). It is worth noting that M92 is of- 



ten considered to be one of the oldest GCs in the Galaxy 
(e.g., Bolte & Hogan 1995). This result gives support to 
the growing notion that the first episodes of star formation 
occurred more or less simultaneously throughout most of 
the Local Group (van den Bergh 1999, 2000; Grebel 2001; 
and references therein). 



On the other hand, recent HST studies of the sparse 
outer-halo GCs Palomar 3, Palomar 4, and Eridanus have 
indicated that they, unlike NGC 2419, might be younger 
than inner-halo GCs of similar metallicity. Since these 
globulars are also commonly considered typical cases of 
loose GCs with predominantly red HBs, they are often 
quoted as "second-parameter" GCs. For instance, Grat- 
ton k, Ortolani (1984) note that the HBs of Pal 3 and Pal 4 
"are too red for [their] metal abundances," thus confirming 
"the need of a second parameter to explain the morphol- 
ogy of the red HBs of outer halo clusters" (see also Stetson 
et al. 1999 and VandenBerg 2000, who explicitly compare 
Pal 3 against M3, on the one hand, and Pal 4/Eridanus 
against M5, on the other hand). 

As widely discussed in the literature, perhaps the pri- 
mary second-parameter candidate is age, though it re- 
mains to be proved that age is the only, or even the most 
important, second parameter (see, e.g., Stetson, Vanden- 
Berg, & Bolte 1996; Sarajedini, Chaboyer, & Demarque 
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TABLE 1 
Observed HB morphology of M3. 
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Fig. 1. — De-reddened CMDs for the internal (panel a) and external (panel b) regions of M3. The region occupied by RR Lyrae variables 
is schematically indicated by boxes. The number of blue (B), variable (V), and red (_R) HB stars is also indicated; note that, for the same 
number of red HB stars, one finds 2.7 times as many blue HB stars in the core than in the cluster outskirts. This sharp difference in HB 
morphology between different radial regions of M3 is real, statistically significant and cannot be caused by any spurious problems with the 
data or the photometry. Two randomly picked CMD simulations for the inner cluster region are shown in the My, (V — I)o plane in panels 
c) and e). Likewise, two randomly picked CMD simulations for the outer cluster region are shown in the My, (B — V)o plane in panels d) 
and f). Note that the total number of stars in the simulations was allowed to randomly fluctuate around the values B + V + R = 174 and 
103 (panels c/e and d/f, respectively), following the Poisson distribution. Note also that the AGB phase is not included in these models. 



1997; Fusi Pecci & Bellazzini 1997; Catelan 2000 and ref- 
erences therein). 

This notwithstanding, that age is the second parame- 
ter has been a frequently adopted assumption in the as- 
tronomical literature, both for Galactic and extragalactic 
GCs (e.g., Sarajedini et al. 1998). Therefore, it is of 
paramount importance to investigate whether the observed 
HB morphology differences are consistent with such a hy- 
pothesis for systems with detected turnoff age differences. 
Arguably, this is also a necessary step before we can be 
in a position to reliably explore the role (if any) of other 
second parameter candidates besides age; after all, if age 
differences can be shown to be compatible with the ob- 
served differences in HB morphology for second-parameter 
pairs, one might invoke Occam's Razor in support of the 
view that, since no other second parameters besides age 
are required, then no other second parameters besides age 
actually exist. 

Original work by Ortolani & Gratton (1989) had previ- 
ously hinted that the "second parameter," in the case of 
Pal 3 and M3, might indeed be age. On the other hand, the 
preliminary HST report by Richer et al. (1996; sec their 
Table 1) supported an age difference of 0.0 Gyr between 
Pal 3 and M3, which — if confirmed — would obviously rule 
out age as the second parameter for this pair. 

More recently, both Stetson et al. (1999) and Vanden- 
Berg (2000) conducted much more detailed comparisons of 
the HST turnoff ages of the outer-halo GCs Pal 3, Pal 4, 
and Eridanus with those of M3 = NGC 5272 (in the case of 
Pal 3) and M5 = NGC 5904 (in the cases of Pal 4 and Eri- 
danus). Curiously, the age differences favored by Stetson 
et al. and by VandcnBcrg differ by ~ 50%, even though 
they are based on the same HST data — with VandcnBerg's 
analysis generally supporting smaller age differences for 
the same second-parameter pairs. (The reasons for such 
differences are discussed in great detail by VandenBerg.) 
Two additional extreme outer-halo GCs, Pal 14 and AM- 
1, were observed by the same team during HST's Cycle 6 
(GO-6512, PI Hesser), but, to the best of our knowledge, 
no results of this investigation have been reported in the 
literature as yet. 

In the previous paper of this series (Catelan 2000), we 
analyzed the outer-halo GCs Pal 4 and Eridanus in the 



context of the "second-parameter" phenomenon. In par- 
ticular, we have quantitatively investigated whether the 
turnoff ages of these outer-halo GCs, as compared to M5, 
are consistent with their HB morphologies (M5 having by 
far the bluer HB). The main conclusion of that study was 
that the HST turnoff age difference, whether as measured 
by Stetson et al. (1999) or by VandenBerg (2000), is too 
small to fully explain the second-parameter effect seen in 
the case of Pal 4/Eridanus vs. M5. 

The purpose of the present paper is to extend Cate- 
lan's (2000) analysis of the outer-halo GCs by carrying 
out a detailed investigation of M3 and Pal 3, which we 
shall assume to be a bona-fide "second-parameter pair" 
(Burbidge & Sandage 1958; Gratton & Ortolani 1984; Or- 
tolani & Gratton 1989; Stetson et al. 1999; VandenBerg 
2000). In particular, our goal is to answer the follow- 
ing question: Assuming that M3 and Pal 3 is a bona-fide 
second-parameter pair, are the HB morphology differences 
between these clusters consistent with the hypothesis that 
age is the second parameter of HB morphology? 

We begin in the next section by describing the HB mor- 
phology of M3 as a function of distance from the cluster 
center, showing that the HB of the cluster is significantly 
bluer in its center than in the outermost regions. The 
CMD of Pal 3, as obtained from the HST-WFPC2 data 
of Stetson et al. (1999), is also presented and discussed; 
previous ground-based photometry of the cluster is also 
briefly addressed. In §3, we build synthetic models of the 
HBs of the two clusters, and present our best-fitting solu- 
tions. In §4, the relative turnoff ages of M3 and Pal 3, as 
determined by Stetson et al. and by VandenBerg (2000), 
are compared with the relative ages that are required to 
account for their different HB types, under the assump- 
tion that age is the second parameter of HB morphology. 
Our concluding remarks, with emphasis on the evidence 
of "local" second parameters leading to the most apparent 
difference in HB morphology between M3 and Pal 3, are 
finally presented in §5. 

2. OBSERVATIONAL DATA 

2.1. HB Morphology of MS 

During the last 10 years we collected and published 
a large photometric data base on M3. Using different 
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Fig. 2. — Illustration of the dependence of several HB morphology parameters on (Mhb) for M3. The mass dispersion parameter was held 
fixed at a value um = 0.023 Mq, appropriate for the inner cluster regions (r < 50"). Each datapoint was obtained from an average of 1000 
HB simulations. Horizontal lines indicate the observed values for M3. In the bottom panel, the solid line indicates B/(B + V + R), whereas 
dotted and dashed gray lines denote V/(B + V + R) and R/(B + V + R), respectively. In the the top panel, BT /{B + V + R) is given, while 
the middle panel shows B0' /B. 



techniques (photographic plates: Buonanno et al. 1994; 
ground-based CCD photometry: Ferraro et al. 1993, Fer- 
raro et al. 1997a, Carrettaet al. 1998; HST-WFPC2: Fer- 
raro et al. 1997b), we secured a huge sample covering the 



entire cluster extension from the very center of the cluster 
up to <~ 7'. On the basis of this dataset we constructed 
one of the largest and most complete luminosity functions 
for evolved stars (brighter than the main-sequence turnoff) 
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Fig. 3. — As in Fig. 2, but for ctm- In these computations, the mean HB mass was held fixed at a value (Mhb) = 0.637 Mq, as found 
appropriate for M3's inner regions (r < 50"). 



ever published for a GC (Rood et al. 1999). 

Using such a dataset we investigated the HB morphology 
in M3 and discovered that it turns out to be significantly 
bluer in the inner region than in the outskirts. While a 
full paper (in preparation) is devoted to accurately dis- 
cussing the HB morphology variation as a function of the 



distance from the cluster center, here we want to point 
out the sharp difference in the HB star distribution be- 
tween the innermost and the outermost regions. Figure 1 
shows the M3 CMDs for the innermost (r < 50", panel a) 
and the outermost (r > 210", panel b) samples. A box 
has been plotted to schematically indicate the location of 
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Fig. 4. — Similar to Fig. 2, but for M3's outer regions (r > 210"); in these computations, the mass dispersion parameter was accordingly 
held fixed at a value um = 0.018 Mq. The upper panel shows the behavior of the Buonanno (1993) parameter, (B2 — R)/(B + V + -R). 



the variable (RR Lyrae-type) stars; the corresponding V 
number is also provided. Figs, lc-f, as well as Figs. 2-5, 
refer to our theoretical modelling of the cluster, and will 
be discussed in §3.1. 

More information on population ratios is summarized 
in Table 1. In column 1, the radial range is given. 
In column 2, the type of data utilized (HST-WFPC2, 
ground-based photographic photometry) is indicated. In 
columns 3, 4, and 5, one finds B, V, and R values, re- 
spectively. The corresponding B : V : R ratios are listed 
in column 6. In columns 7 through 9, the values of the 
Mironov (1972), Lee-Zinn (Zinn 1986; Lee et al. 1994), 
and Buonanno (1993) HB morphology parameters, respec- 
tively, are given. The last three columns provide the values 
of some additional HB morphology parameters defined by 
Catelan et al. (2001); we briefly recall that B2' represents 
the number of HB stars redder than (V — I)o = —0.02, 
whereas BT and B0' represent the number of blue HB 



stars redder than (V — I)o = +0.07 and 0, respectively. 

As can be seen from Fig. 1, in selecting different groups 
along the HB we followed the criteria described in Buo- 
nanno et al. (1994) and Ferraro et al. (1997a): in partic- 
ular, both "Extreme Blue" (EB), with V > 17, and "Ex- 
treme Red" (ER), with [B-V) > 0.6 or (V - I) > 0.8, 
HB stars have been excluded from the B and R number 
counts, respectively. In fact, these stars might be "non 
genuine" HB stars, in the sense that they may derive from 
different (non-canonical) phenomena: ER objects might 
represent the progeny of the large blue straggler popula- 
tion discovered in this cluster (Ferraro et al. 1993, 1997b); 
and EB stars could derive from (say) late collisonal merg- 
ers (Bailyn & Iben 1989). In any event, ER and EB stars 
represent a very minor fraction of the total HB population 
of M3 — 3.3% and 2.6%, respectively — so that the results 
of this paper are not seriously affected by our not having 
included them in the number counts reported in Table 1. 
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Fig. 5. — As in Fig. 4, but for <tm- In these computations, the mean HB mass was held fixed at a value (Mhb) = 0.646 Mq, as found 
appropriate for M3's outer regions (r > 210"). 



As can be seen from the data plotted in Figs. la,b and 
summarized in Table 1, the detected difference in the M3 
HB morphology is sharp and significant, ft is difficult to 
attribute it to any spurious effects (due to crowding, etc.) 
or to any additional source of error associated with the 
photometry. Furthermore, the limiting magnitude of the 
M3 photometry is several magnitudes below the instability 
strip level in Figs. la,b, ruling out significant differential 
incompleteness in HB star detection from one side of the 
HB to the other. In our view this represents clear-cut 
evidence that a mechanism able to modify the star distri- 
bution along the HB is at work in M3 and is more effective 
in the central region of the cluster. The problem is surely 
worthy of further investigation since anomalous radial be- 
haviour has been seen before in other stellar populations 
in this cluster (see the case of the blue straggler radial dis- 
tribution discussed in Ferraro et al. 1993 and Ferraro et 
al. 1997b). 



2.2. HB Morphology of Palomar 3 

The dereddened V, {V-I) Pal 3 CMD, based on HST- 
WFPC2 data kindly provided to us by Dr. P. B. Stetson, 
is shown in Figs. 6a, b. RR Lyrae variables are indicated 
by circled dots. The adopted reddening comes from the 
Schlegel, Finkbeiner, & Davis (1998) dust maps, which 
give E{B — V) = 0.042 mag at the position of Pal 3, and 
the assumption that E(V - I)/E{B-V) ~ 1.3 (see Stet- 
son et al. 1999 and references therein). The adopted red- 
dening is in excellent agreement with the value tabulated 
by Harris (1996) in the June 1999 update of his catalog, 
E(B — V) = 0.04 mag — an estimate which is completely 
independent of the Schlegel et al. maps. For reference pur- 
poses, a vertical dotted line representing the mean dered- 
dened color of the combined HBs of Pal 4 and Eridanus 
[(V — I)o ~ 0.78 mag; see Catelan 2000] is overplotted on 
the Pal 3 CMD of Figs. 6a,b. In Fig. 6b, the zero-age HB 
for Z = 0.001 is overplotted on the data, as is an evolution- 
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ary track for a 0.656 Mq star, after appropriately applying 
a shift in V by 19.67 mag. 

Mean colors and magnitudes for the Pal 3 RR Lyrae 
variables were obtained by constructing their light curves. 
This was achieved by matching the Borissova et al. (1998) 
and Stetson et al. (1999) photometries, which allowed 
light curve types, pulsation periods and mean properties 
to be evaluated (J. Borissova 2000, priv. comm.). While 
these cannot yet be considered the definitive light curves 
for the Pal 3 variables, they certainly represent a substan- 
tial improvement in comparison with Stetson et al. , where 
only a global "mean" color/magnitude was provided for the 
cluster variables. Importantly, it seems particularly clear 
that: i) All of the Pal 3 RR Lyrae variables pulsate in 
the fundamental (ab) mode; ii) All of the Pal 3 RR Lyrae 
variables are redder than (V — I)o ~ 0.41 mag: this is 
the blue limit for RRab pulsation in M3, the prototypical 
Oostcrhoff type I cluster (Bakos & Jurcsik 2000; Corwin & 
Carney 2001). Note that RRab's are even redder in Oost- 
erhoff type II globulars (e.g., Bono, Caputo, & Marconi 
1995). 

From Figs. 6a, b, it is immediately obvious that the HB 
of Pal 3 is much bluer than the HBs of Pal 4 and Eri- 
danus (whose mean combined color is indicated, in these 
plots, by a vertical dotted line, from Catelan 2000). How- 
ever, blue-HB stars appear to be entirely absent from all 
these outer-halo globulars. Pal 3 appears to be just slightly 
more metal-poor than Pal 4 and Eridanus (see Stetson et 
al. 1999 for a detailed discussion). The small metallicity 
difference between Pal 3 and Pal 4/Eridanus suggests that 
metallicity differences alone cannot explain the difference 
in HB color between these clusters as seen in Figs. 6a, b. 
As a matter of fact, Burbidge & Sandage (1958), in their 
analysis of Pal 3 and Pal 4 (which they call the "10 h " 
and the "ll h " clusters in their paper, respectively) had 
already pointed out that the Pal 3 CMD is "more nor- 
mal" than that of Pal 4 "because the HB extends to bluer 
colors." In spite of that remark, and the existence of sev- 
eral RR Lyrae variables in the cluster (unlike in Pal 4 and 
Eridanus) notwithstanding, Pal 3 has usually been consid- 
ered as good an example of a "second-parameter" cluster 
as any of the other rcd-HB globulars in the outer halo. 
Our qualitative analysis of the CMD presented in Fig. 6 
suggests, in line with Burbidge & Sandage, that Pal 3 is 
a much less extreme second-parameter cluster than Pal 4 
and Eridanus. As we shall see below, a detailed quantita- 
tive analysis confirms that this is indeed the case. 

As already pointed out, the RR Lyrae population of 
Pal 3 has not been extensively investigated to date. How- 
ever, the currently known RR Lyrae variables in the clus- 
ter all appear to be of the ab-type — i.e., fundamental pul- 
sators (Borissova et al. 1998; Borissova 2000, priv. comm.; 
Stetson et al. 1999 and references therein). RRab Lyrae 
variables are well known to be redder than RRc variables, 
implying that the HB of Pal 3, which lacks either blue HB 
or RRc stars, has a narrow dispersion in color. Note that 
this argument independs on the reliability of the individ- 
ual RR Lyrae colors. The mean color of the Pal 3 HB 
stars, variable and non-variable, appears to coincide with 
the red edge of the instability strip, at (V — I)o « 0.58 mag 
(see Figs. 6a, b). We shall argue, in the next section, that 
the three fairly bright stars at the red extension of the 



HB, at (V — I)q w 0.8 — 0.85 mag, are red HB stars on 
their way to the asymptotic giant branch (AGB), even 
though the average expected number of such evolved stars 
is ~ 1. Taking this into consideration in the number 
counts, we conclude that the most representative num- 
ber ratios and mean color describing the HB morphol- 
ogy of Pal 3 are B : V : R = 0.00 : 0.50 : 0.50 and 
({V - I)o)hb = 0.583 ± 0.028 mag, respectively. 

Gratton & Ortolani (1984) suggest that "one or two 
stars may be present on the blue side of the RR Lyrae 
strip" of Pal 3. Even though these two stars are very likely 
cluster members (K. Cudworth 2000, priv. comm.), inde- 
pendent photometry by K. Cudworth and by S. Ortolani 
(2000, priv. comm.) shows that they are both variable 
stars (at least one of which almost certainly being an RR 
Lyrae variable), and not bluc-HB stars as had originally 
been suggested. In fact, Ortolani has imaged a Pal 3 field 
of about 5' x 5' using the ESO-NTT in 1993, not finding 
any blue HB star in the cluster. Ortolani points out that 
Pal 3, in his images, seems to have a diameter much smaller 
than his observed field — namely, no more than ~ 1.5'. In- 
dependently, Borissova (2000, priv. comm.) has reduced 
B and V images taken with the Rozhen telescope and cov- 
ering a field of 5.6' x 5.6' centered on the Pal 3 center. The 
limiting magnitude, in her study, is 23 mag in V (compare 
with Fig. 6). Her analysis of 25 images taken in B and 
20 images taken in V also confirms that Pal 3 contains 
no blue HB stars. Given that the Pal 3 half-light radius is 
only rh ~ 40" (Harris 1996), the conclusion that Pal 3 does 
entirely lack blue HB stars seems inescapable. Hence the 
absence of blue HB stars in the HST CMD of Figs. 6a,b 
is not an artifact of the relatively small field of view of 
WFPC2 (compared to the ground-based ones). 

Of course, it would obviously be desirable to have a 
larger sample of HB stars. Note however that even pho- 
tometry covering the entire cluster out to its tidal radius 
is not likely to yield a large sample. Pal 3's tidal radius is 
r t = 4.81' (Harris 1996); a simple integration assuming a 
King profile suggests that the number of HB stars in the 
range not covered by the Borissova (2000, priv. comm.) 
or Ortolani (2000, priv. comm.) photometries mentioned 
in the previous paragraph — i.e., in between about 3' and 
5' — is very small, of order a few percent of the total inside 
r ~ 3' only, or less than one HB star. On the other hand, it 
is possible to estimate the "total" number of Pal 3 HB stars 
from Eq. (14) in Rcnzini & Buzzoni (1986), which gives the 
expected number of stars in any given post-main sequence 
stage as a function of the duration of that stage and of the 
total luminosity of the sampled population. Adopting a 
lifetime of 10 8 yr for the HB, and a total integrated mag- 
nitude My = —5.6 mag for Pal 3 (from Harris 1996), we 
find that f ot w 30. This number seems surprisingly 
large in comparison with the number of HB stars detected 
in the cluster so far. We have no explanation at present 
for this discrepancy, but believe that it is highly unlikely 
that the answer will be found in the form of blue HB stars 
outside the currently observed cluster regions. 

3. THEORETICAL FRAMEWORK 

In this section, our goal is to determine the mean HB 
mass and HB mass dispersion parameters ((Mhb), om ) in 
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Table 2. HB morphology of M3: simulation results. 



Range 


(B) : (V) : (R) 


(—) 


( B ~ R ) 


/ B2-R \ 
\B+V+RI 


( B2 ' ) 
\b+v+ri 


\b+v+ri 


(f) 


r < 50" 


0.472 : 0.334 : 0.194 


0.709 


0.279 


-0.052 


0.917 


0.238 


0.755 




(0.039) (0.037) (0.030) 


(0.043) 


(0.059) 


(0.044) 


(0.021) 


(0.032) 


(0.047) 


r > 210" 


0.298 : 0.443 : 0.259 


0.536 


0.039 


-0.218 


0.984 


0.209 


0.906 




(0.044) (0.047) (0.043) 


(0.066) 


(0.073) 


(0.051) 


(0.012) 


(0.041) 


(0.052) 



the simulations that best account for the observed char- 
acteristics of the M3 (§3.1) and Pal 3 (§3.2) CMDs in the 
HB region. 

The basic theoretical framework adopted (HB and RGB 
evolutionary tracks, color transformations) is the same as 
in Catelan (2000 and references therein), and need not be 
repeated here; we refer the interested reader to that pa- 
per for further information. Our employed HB synthesis 
code, SINTDELPHI, is the latest version of the Catelan et al. 
(1998 and references therein) code, which includes several 
numerical improvements and has now been translated to 

DELPHI. 

Our most important working hypothesis is that age is 
the second parameter of HB morphology. We also assume, 
in line with previous authors, that M3 and Pal 3 repre- 
sent a "bona fide" second-parameter pair, thus adopting 
the same metallicity for the two clusters — which appears 
to be supported by the available data (see Stetson et al. 
1999 and references therein). In particular, Armandroff, 
Da Costa, & Zinn (1992) have measured a metal abun- 
dance for Pal 3 using the Can triplet method, and found 
[Fe/H] = -1.57 ± 0.19 dcx on the Zinn & West (1984) 
scale. This value is to be compared with [Fe/H] = —1.57 
for M3, as provided by Harris (1996), on the same scale. 
It should be noted that the implied identity between the 
Pal 3 and M3 metal content will eventually have to be 
placed on a firmer footing with the use of high-dispersion 
spectroscopy; in particular, a much smaller error bar than 
provided for Pal 3 in the Armandroff et al. study will be 
needed to verify the assumption that Pal 3 and M3 have 
indeed the same metallicity. 

Again as in Catelan (2000), we assume that the under- 
lying zero-age HB mass distributions can be characterized 
by Gaussian deviates (see Catelan et al. 1998 for a detailed 
discussion). 

3.1. The Case of MS 

Extensive grids of synthetic HBs have been computed 
aiming at estimating the optimum parameters (Mhb) 
(mean mass) and om (mass dispersion) required to re- 
produce the observed HB morphology of M3 at both its 
bluest and reddest extremes (see §2.1). 

Random "observational scatter" has also been included, 
according to the following exponential "laws" that were 
fitted to the observational data: 

<j b = 0.004 exp(M B - 3.7) + 0.025, 
cry = 0.004 exp(M v - 3.8) + 0.018, 



a i = 0.004 exp(M/ - 3.8) + 0.018. 

We find that these formulae reproduce the errors in the 
photometry over the whole range of magnitudes more ac- 
curately than do similar ones adopting Robertson (1974) 
profiles. 

The width of the theoretical instability strip was as- 
sumed to be A log T e ff = 0.085, which leads to a very good 
agreement with the empirical blue and red edges of the M3 
RR Lyrae strip as reported by Sandage (1990), Bakos & 
Jurcsik (2000), and Corwin & Carney (2001). Moreover, 
the adopted width also leads to a nice agreement between 
the model predictions and the observed red edge of the 
instability strip in the case of Pal 3, which is particularly 
important for our present purposes (see below). 

We have computed series of 1000 synthetic HBs assum- 
ing an overall number of HB stars B + V + R = 174 and 
103 for the "blue" and "red" HB cases, respectively. These 
numbers correspond to the actual number counts in the in- 
ner and outer regions of the cluster, respectively (see Ta- 
ble 1), and were allowed to fluctuate following the Poisson 
distribution. 

The number ratios B : V : R provided in Table 1, as 
well as the dispersion in color in the observed HBs, allowed 
us to iteratively fine-tune the model simulation parameters 
{(Mhb), cm) until the best combinations were obtained. 
Figures 2 and 3 illustrate the dependence of the employed 
HB morphology parameters on the mean HB mass and HB 
mass dispersion, respectively, for the inner regions of the 
cluster. Figures 4 and 5 are analogous and refer to the 
outer regions of M3. Such figures were also used to obtain 
approximate estimates of the corresponding errors in the 
best-fitting values. We found that the following provides 
the best match to the observed HB morphology of M3 in 
its inner regions: 

(M HB ) = 0.6370 ± 0.0025 M Q , a M = 0.023 ± 0.004 M Q ; 

whereas, for the outer regions, the best match is found for: 

(M HB ) = 0.6457 ± 0.0025 M , a M = 0.018 ± 0.003 M Q . 

Such combinations lead, in the mean, to the HB morphol- 
ogy parameters displayed in Table 2 (where the numbers in 
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Fig. 6. — Panel a) shows the de-reddened HST-WFPC2 CMD for Palomar 3. As indicated, a reddening of E(V — I) = 0.055 mag has been 
assumed (based on Schlcgel et al. 1998). RR Lyrae variables are indicated by circled dots. The vertical dotted line indicates the mean color 
of the bulk of the HB populations in the outer-halo GCs Pal 4 and Eridanus, ((V — I)o) = 0.78 mag (see Catelan 2000). Panel b) is identical 
to panel a), except that in this case a zero-age HB for Z = 0.001 (dashed line) is ovcrplotted on the data, as is an evolutionary track for an 
HB star with the same metallicity and M = 0.656 Mq, after a suitable shift in V by 19.67 mag. In panels c) through h), a random sample of 
CMD simulations for this cluster ("best-fitting case") is provided for comparison. Note that the simulations do not include the AGB phase. 



parentheses represent the standard deviation of the mean 
over the corresponding set of 1000 simulations). Note the 
nice agreement between the observed (Table 1) and theo- 
retical (Table 2) number counts. We also emphasize the 
agreement between observed and simulated values of the 
Mironov (1972) parameter which, unlike the Lee-Zinn 
parameter, is immune to possible extant uncertainties in 
the RR Lyrae number counts in the innermost regions of 
M3 observed only with HST. Note also that the best-fitting 
mass dispersion appears to be slightly larger in the inner 
than in the outer regions of M3, although it will imme- 
diately be recognized that this result does not bear great 
statistical significance. 

Four simulations were randomly drawn from the two sets 
of 1000 simulations whose mean HB morphology param- 
eters best reproduce the observed M3 ones over the two 
radial regions considered. The corresponding synthetic 
CMDs are plotted in Figs. lc,e [My, (V — I)o plane] and 
Figs. ld,f [My, (B — V)o plane] for the inner and outer re- 
gions, respectively. For each simulation, the total number 
of HB stars, as well as the relative fractions B : V : R, 
are indicated. The simulations do not include the AGB 
phase. 

Our simulations allow us to address the question 
whether the difference between the "inner cluster" and 
the "outer cluster" HB populations is significant or not. 
We have proceeded as follows. For the set of simulations 
which best reproduces, in the mean, the observed number 
counts in the inner regions of M3, we have selected those 
for which V > B (i.e., more RR Lyrae variables than blue- 
HB stars) — a characteristic of the outer regions of M3. 
We find that only about 3% of the "inner cluster" syn- 
thetic HBs have "outer cluster" characteristics, with more 
RR Lyrae variables than blue-HB stars. Likewise, about 
6% of the "outer cluster" simulations have "inner clus- 
ter" characteristics, in the sense that they show B > V. 
Moreover, since the number of RR Lyrae and (especially) 
red HB stars are essentially the same for the inner and 
outer samples (cf. Table 1 and Figs. la,b), the difference 
arises mainly from the number of blue HB stars. If the two 
populations are indeed the same our best estimate of the 
number of blue HB stars in each sample is 55.5. If only 
Poisson statistics were acting the observed inner sample 
is a 3.5-0- deviation upward occuring simultaneously with 
a 3.5-ct deviation downward in the outer sample. Other 
scenarios are similarly improbable on a purely statistical 
basis. Hence the difference in HB morphology between the 
inner and outermost regions of M3 is clearly of very high 
statistical significance. 

3.2. The Case of Pal 3 

The HB of Pal 3 that we set out to simulate theoreti- 
cally, from the HST observations of Stetson et al. (1999), 
completely lacks blue HB stars (see Figs. 6a, b). Inspec- 
tion of Figs. 6a,b clearly shows that the bulk of the HB 



stars cluster around the red edge of the instability strip, 
which is located at a color (V — I)o — 0.58 mag. With 
the exception of the three bright, likely evolved stars at 
(V — 7)o — 0.8 — 0.85 mag, the dispersion in color among 
the HB stars appears to be surprisingly small, especially 
given the strong sensitivity of temperature to mass varia- 
tions at these intermediate colors (e.g., Sweigart, Renzini, 
& Tornambc 1987). This is confirmed by the fact that 
all known or suspected RR Lyrae variables in the cluster 
are fundamental (i.e., ab-type) pulsators (Borissova et al. 
1998; Stetson et al. 1999; Borissova 2000, priv. comm.). 

The adopted procedure in the synthetic computations 
is entirely analogous to the one already described in the 
case of M3. "Observational scatter" has been randomly in- 
cluded following exponential "laws" that were fitted to the 
HST data. The best-fitting error formulae are as follows: 



cry = 0.0008 exp(My - 2.5) + 0.0085, 



07 = 0.0007 exp(Mj - 1.5) + 0.0075. 



In these simulations, the number of HB stars was al- 
lowed to fluctuate around the observed number, B + V + 
B = 16, following the Poisson distribution. 

As a result, we find that the following parameters lead, 
in the mean, to the best overall match between the ob- 
served and simulated Pal 3 HBs: 



(M HB > = 0.661 ± 0.004 M Q , a M = 0.0035±8;gg|g M< 



The mean number ratios and (V — I)o color for this set of 
simulations is given in Table 3. A sample of six simula- 
tions, drawn randomly from the pool of 1000 models, is dis- 
played in Figs. 6c-h. Comparison between Figs. 6a, b and 
Figs. 6c-h confirms the good agreement between the sim- 
ulations and the observed CMD. Note, in particular, that 
the mean computed color is in excellent agreement with 
the one based on the HST data (§2.2). The only aspect of 
the observed CMD which may perhaps be argued not to 
be well matched by the simulations is the slope of the HB. 
However, much of the sloping nature of the Pal 3 HB is due 
to a single red HB star, at around (V — J)o ~ 0.65 mag, 
that appears to be significantly fainter than the remaining 
red HB stars; this is certainly not expected on the basis of 
our simulations (see also VandenBerg 2000). The compar- 
ison between the models and the observations carried out 
in Fig. 6b suggests that, with the exception of this intrigu- 
ingly faint star, there is no obvious discrepancy between 
the evolutionary models and the observations. 
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Fig. 7. — Illustration of the dependence of several HB morphology parameters on (Mhb) for Pal 3. The mass dispersion parameter was held 
fixed at a value cm = 0.0035 Mg, as found appropriate for the cluster. Each datapoint was obtained from an average of 1000 HB simulations. 
Horizontal lines indicate the observed values. In the bottom panel, the solid line indicates B/(B + V + R), whereas dotted and dashed gray 
lines (which actually overlap) denote V/(B + V + R) and R/(B + V + R), respectively. In the top panel, the standard deviation in color of 
the HB distribution, a(V — /)hb, is given; the middle panel shows the mean HB color, {(V — -Oo.hb)- 



Table 3. HB morphology of Palomar 3 



(B) : (V) : (R) (- L 



K ) «V-/)o,iib> o(V-I)» 



Observed 0.00 : 0.50 : 0.50 -0.50 
(0.00) (0.23) (0.2.3) (0.2.3) 



Simulations 0.00 : 0.49 : 0.51 

/n nn \ /n ill (n 191 



-0.51 

I Pi 1 Q \ 
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Fig. 8. — As in Fig. 7, but for cm- In these computations, the mean HB mass was held fixed at a value (Mhb) = 0.661 Mq, as found 
appropriate for Pal 3. 

Fig. 9. — Same as in Fig. 6, except that in panels c) through h) a random sample of CMD simulations for Pal 3 is provided in the wide, 
M3-like (cm = 0.02 Mq) mass dispersion case. Note that, unlike in the observed CMD (panels a and b), all of the simulations displayed here 
contain one or more stars bluer than (V — I)o = 0.40 mag, including at least one bona-fide blue-HB star. 



Figures 7 and 8 show the run of several HB morphology HB number counts do play an important role when choos- 
parameters for Pal 3 as a function of the adopted (Mhb) ing the best-fitting (Mhb) value (Fig. 7, bottom panel), 
and gm, respectively. It will be noted that, whereas the the mean HB color, ((V — /)o,hb), is (not unexpectedly) 
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Fig. 11. — Fraction of synthetic "Palomar 3" HB simulations containing at least one star bluer than (V — I)o = 0.4 mag (open squares) or 
(V — I)o = 0.45 mag (filled circles). In these simulations, the mean mass was held fixed at the best-fitting value for Pal 3, (Mhb) = 0.661 M@, 
but the mass dispersion ctm was allowed to vary over a wide range. The acceptable M3 mass dispersion range (see §3.1) is schematically 
shown; fr M could actually be larger than indicated for M3 if any of the rejected "ER" or (especially) "EB" stars (§2.1) turn out to be "bona 
fide" HB stars. This figure provides an argument, based solely on its lack of sufficiently blue stars, that the Pal 3 mass dispersion on the HB 
is smaller than M3's at the > 2-cr level. Note that each data point corresponds to a scries of 10,000 Monte Carlo simulations. 



another extremely important quantity constraining the 
(Mhb) solution. On the other hand, Fig. 8 shows that the 
traditionally employed number counts B, V, R do not pose 
useful constraints on <7m, and neither does the mean HB 
color; in order to estimate the best-fitting ctm value, the 
standard deviation in the HB color, a(V — I)hb, is clearly 
the most important quantity (Fig. 8, upper panel). It is 
important to note, at this point, that in order to maximize 
the sensitivity of a(V — I)hb to the HB mass dispersion, 
evolved stars were avoided both in the models and in the 
observations. In the latter case, the three evolved Pal 3 
HB stars were not included when computing a(V — /)hb; 
in the former, we have decided to simply exclude from the 
computations any star approaching helium exhaustion in 
the core, or (equivalently) the final 10% of HB evolution. 
For these reasons, it is not surprising that a(V — I)hb 
provides a very good diagnostic of mass dispersion on the 
zero- age HB. 

The displayed models also shed light on the nature of the 
three red and bright stars lying to the red of the "main" 
Pal 3 HB distribution. For instance, one of the displayed 
simulations, labeled #942 (Fig. 6c), shows two bright and 
red HB stars, bearing great resemblance to the observed 
distribution. These stars are approaching their final stages 



of HB evolution, on their way up to the AGB. As shown in 
the figure, the existence of such bright and red HB stars is 
not unexpected from our models (see also Catelan 2000), 
although we find that the typical number to be expected 
in a sample the size of Pal 3's (16 HB stars) is ~ 1. This 
is due to the relatively short duration of the final stages 
of HB evolution close to core helium exhaustion. Though 
suggestive, the three observed bright and red HB stars in 
Pal 3 cannot be considered a serious discrepancy between 
theory and observations in their own right. 

One particularly interesting result of our model compu- 
tations is the extremely small value of <jm that we have ob- 
tained for Pal 3. Such a small mass dispersion is largely a 
consequence of: i) The lack of blue HB stars in the cluster; 
ii) The lack of first-overtone RR Lyrae variables (RRc's). 
One may, however, question the statistical significance of 
this result, given that the Pal 3 field that we are modelling 
has a mere 16 HB stars. For this reason, we have computed 
a new set of simulations for Pal 3 with the same value of 
(Mhb) = 0.661 M Q as derived above, but changing the 
mass dispersion to a value that appears more appropriate 
for M3, namely: a M = 0.02 M (see §3.1). 

A random sample of six simulations drawn from the so 
computed pool of 1000 models is displayed in Fig. 9. It 
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is readily apparent that the typical Pal 3 simulation with 
an M3-like mass dispersion on the HB shows an HB distri- 
bution which does not provide a good match to the Pal 3 
CMD. In particular, we find that only 12.3% of these sim- 
ulations entirely lack blue-HB stars, as in the HST CMD. 
Even in these cases, however, good morphological agree- 
ment with the observed Pal 3 CMD is generally not ob- 
tained. To illustrate this, a random sample of six simula- 
tions drawn from the (already small) subset of simulations 
with B — is shown in Figs. lOc-h. None of the examples 
shown has a narrow HB color dispersion resembling the 
observed one (compare with Figs. 10a, b), and RRc stars 
are often present in these simulations (as indicated by the 
instability strip being populated to bluer colors than in 
Figs. 6c-h, where very few RRc's are present). As a mat- 
ter of fact, the bluest of the stars plotted in Figs. 10a, b 
has a color (V — I)o = 0.46 mag, whereas all of the sim- 
ulations displayed in Figs. lOc-h contain stars bluer than 
(V — I)q = 0.40 mag. From the entire pool of simulations 
for this M3-like mass dispersion case, we indeed find that 
as many as 98.4% contain one or more stars bluer than 
(V — I)o = 0.40 mag; this percentage increases to 99.1% 
if one considers (V — I)q = 0.45 mag as the "cutoff color" 
on the blue side. We conclude that the small HB mass 
dispersion that we found for Pal 3 is most likely an in- 
trinsic characteristic of Pal 3, and not a chance result due 
to statistical fluctuations. This conclusion is also strongly 
supported by the upper panel of Fig. 8, which clearly shows 
not only that the HB mass dispersion of Pal 3 is smaller 
than 0.01 M Q — but also that it is, perhaps somewhat sur- 
prisingly, entirely consistent with zero. Such a small mass 
dispersion for the Pal 3 HB is due to the fact that the 
HB evolutionary tracks, at the Pal 3 mean color, already 
encompass a significant range in color at this metallicity. 
This is clearly demonstrated by Fig. 6b, which shows an 
evolutionary track from Catelan et al. (1998) for an HB 
star with a mass very similar to Pal 3's best-fitting (Mhb) 
overplotted on the observed CMD. 

To put these results on a firmer statistical basis, we have 
computed many additional sets of simulations in which the 
mean HB mass was held fixed at the "optimum" value for 
Pal 3, namely (Mhb) = 0.661 M Q , but allowing the mass 
dispersion um to vary between 0.0 M Q ("single HB evo- 
lutionary track case") and 0.035 M . For each um value, 
10,000 Monte Carlo simulations with 16 HB stars in each 
(as observed), plus Poisson noise, were computed, and the 
number of those simulations showing one or more HB star 
bluer than the "cutoff colors" (V — I)o = 0.40 mag and 
(V — I)o = 0.45 mag obtained. The results of these sim- 
ulations arc summarized in Fig. 11, where the percentage 
of simulations containing one or more star bluer than the 
indicated cutoffs is plotted as a function of the mass dis- 
persion o~m ■ The range of mass dispersion values found for 
M3 (from §3.1) is also schematically indicated; note that 
these may be underestimates, in case any of the rejected 
"ER" and "EB" stars (§2.1), particularly the latter, turn 
out to be "bona fide" HB stars. 

Figure 11 clearly confirms that the lack of bluer HB 
stars in the Pal 3 HST-WFPC2 CMD implies that the 
HB mass dispersion in this cluster is smaller than M3's 
at the > 2-cr level. This is true even when the cutoff 
color is placed at (V — I)o = 0.40 mag, which is 0.06 mag 



bluer than the reddest point on the HST-WFPC2 CMD. 
The choice of such a significantly bluer cutoff color was 
motivated by the fact that, in spite of recent efforts, the 
Pal 3 RR Lyrae light curves are still not very well de- 
termined, so that their individual colors remain somewhat 
uncertain. However, we again emphasize that all RR Lyrae 
stars found to date in Pal 3 appear to be fundamental pul- 
sators. From Fig. 2 in Bakos & Jurcsik (2000) and from 
Corwin & Carney (2001) it is clear, in turn, that all RRab 
stars in M3 (a prototypical Oostcrhoff type I cluster) are 
redder than (V — I)o ~ 0.41 mag; RRab's are even redder 
in Oosterhoff type II globulars (e.g., Bono et al. 1995). 

In conclusion, it appears clear that the mass dispersion 
on the Pal 3 HB is very small, being not only intrinsically 
smaller than M3's but also indistinguishable from zero. 

4. ESTIMATING RELATIVE AGES FROM DIFFERENCES IN 
THE HB MORPHOLOGY 

In order to estimate the relative ages required to pro- 
duce the relative HB types of M3 and Pal 3, we follow the 
same approach described in Catelan (2000) in the case of 
M5 and Pal 4/Eridanus. As in that paper, we shall evalu- 
ate the effects of an age-dependent mass loss on the RGB, 
as implied by the several different analytical formulae dis- 
cussed in Catelan's Appendix. As argued by Catelan, at 
the present time it does not appear possible to strongly 
argue in favor of any of his recalibrated equations over the 
others, so that a safer approach is to use all of them simul- 
taneously whenever an estimate of the amount of mass loss 
on the RGB is needed. We briefly recall here that equa- 
tions (Al) through (A4) of Catelan are "generalized," em- 
pirically recalibrated variations of the analytical mass loss 
formulae previously suggested by Reimers (1975a, 1975b), 
Mullan (1978), Goldberg (1979), and Judge & Stencel 
(1991), respectively. In this paper, RGB mass loss was 
estimated on the basis of the RGB models of VandenBerg 
et al. (2000) for a chemical composition [Fe/H] = —1.54, 
[a/Fe] = +0.3. 

We follow the same approach described in Catelan 
(2000) to measure the age difference between the two clus- 
ters, as required from their HB morphologies under the 
assumption that age is the second parameter. Given the 
variation in HB morphology as a function of radial po- 
sition in the case of M3, we are forced to study the M3 
inner regions (bluer HB) and outer regions (redder HB) as 
limiting cases for deriving the HB morphology-based age 
difference between M3 and Pal 3; we assume that the most 
appropriate age difference estimate will lie somewhere in 
between these limits. In particular, the "integrated" HB 
morphology of M3 can be characterized by the number 
ratios B : V : R = 0.39 : 0.40 : 0.21, which would corre- 
spond to a synthetic HB falling in between the two limiting 
cases considered above (compare with Tables 2 and 3). 

For each of the Catelan (2000) mass loss formulae, and 
also for Reimers' (1975a, 1975b) formula, the resulting r\ 
values and "HB morphology ages" for M3 and Pal 3 are 
given in Table 4 ( "red HB" case for M3) and Table 5 ( "blue 
HB" case for M3). These results are also shown in Fig. 12, 
where one can find our theoretical age difference vs. ab- 
solute age curves for each mass loss formula (see labels) 
both in the "red HB" and "blue HB" limits for M3 (upper 
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TABLE 4 

HB MORPHOLOGY AGES FOR M3 ( "RED" HB CASE) AND PAL 3. 
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0.74 


3.165 


16.24 


0.76 


15.0 


0.263 


14.30 


0.70 


1.656 


14.40 


0.60 


4.337 


14.37 


0.63 


13.0 


0.349 


12.41 


0.59 


2.291 


12.47 


0.53 


5.930 


12.46 


0.54 


11.0 


0.462 


10.54 


0.46 


3.133 


10.59 


0.41 


8.057 


10.58 


0.42 


9.0 


0.636 


8.73 


0.27 


4.546 


8.78 


0.22 


11.561 


8.77 


0.23 




Equation (A3) 




Equation (A4) 




Equation (A5) 




21.0 


4.559 


19.89 


1.11 


0.644 


20.00 


1.00 


1.287 


19.95 


1.05 


19.0 


6.208 


18.00 


1.00 


0.914 


18.13 


0.87 


1.790 


18.07 


0.93 


17.0 


8.249 


16.14 


0.86 


1.277 


16.28 


0.72 


2.440 


16.22 


0.78 


15.0 


10.885 


14.28 


0.72 


1.792 


14.42 


0.58 


3.316 


14.35 


0.65 


13.0 


14.250 


12.37 


0.63 


2.518 


12.50 


0.50 


4.479 


12.43 


0.57 


11.0 


18.506 


10.52 


0.48 


3.520 


10.62 


0.38 


5.986 


10.56 


0.44 


9.0 


25.219 


8.74 


0.26 


5.210 


8.78 


0.22 


8.456 


8.76 


0.24 



TABLE 5 

HB MORPHOLOGY AGES FOR M3 ("BLUE" HB CASE) AND PAL 3. 



tg 1) tg Atg 1) tg Atg 1) tg Atg 



(M3) 




(Pal 3) 






(Pal 3) 






(Pal 3) 








Reimers 




Equation (Al) 




Equation (A2) 




21.0 


0.120 


19.30 


1.70 


0.676 


19.44 


1.56 


1.830 


19.40 


1.60 


19.0 


0.161 


17.49 


1.51 


0.935 


17.65 


1.35 


2.507 


17.61 


1.39 


17.0 


0.212 


15.72 


1.28 


1.278 


15.87 


1.13 


3.391 


15.83 


1.17 


15.0 


0.279 


13.92 


1.08 


1.754 


14.06 


0.94 


4.594 


14.03 


0.97 


13.0 


0.366 


12.09 


0.91 


2.403 


12.18 


0.82 


6.221 


12.16 


0.84 


11.0 


0.481 


10.29 


0.71 


3.258 


10.37 


0.63 


8.379 


10.36 


0.64 


9.0 


0.657 


8.58 


0.42 


4.694 


8.65 


0.35 


11.938 


8.64 


0.36 




Equation (A3) 




Equation (A4) 




Equation (A5) 




21.0 


5.053 


19.28 


1.72 


0.714 


19.46 


1.54 


1.426 


19.38 


1.62 


19.0 


6.746 


17.47 


1.53 


0.993 


17.68 


1.32 


1.945 


17.58 


1.42 


17.0 


8.838 


15.68 


1.32 


1.368 


15.90 


1.10 


2.614 


15.80 


1.20 


15.0 


11.531 


13.88 


1.12 


1.898 


14.10 


0.90 


3.513 


13.99 


1.01 


13.0 


14.948 


12.02 


0.98 


2.642 


12.23 


0.77 


4.699 


12.12 


0.88 


11.0 


19.246 


10.26 


0.74 


3.661 


10.41 


0.59 


6.226 


10.32 


0.68 


9.0 


26.040 


8.59 


0.41 


5.380 


8.66 


0.34 


8.732 


8.63 


0.37 



panel and lower panel, respectively). 

VandenBerg (2000), on the basis of deep WFPC2 pho- 
tometry for Pal 3, has recently found a turnoff (TO) age 
"of ss 11.5 Gyr for Pal 3, which is to be compared with 
12 - 12.5 Gyr for M3" (see his §2.4.2)— implying a favored 
age difference for the M3-Pal 3 pair of « 0.5 — 1 Gyr. Com- 
paring such a TO age difference with the HB morphology- 
based results displayed in Fig. 12, we find that the relative 
TO ages favored by VandenBerg are largely consistent with 
the relative HB types of the M3 and Pal 3 HBs. 

On the other hand, Stetson et al. (1999) state that "the 
[deep WFPC2] data are consistent with an age difference 
of some 2 Gyr between M3 and Pal 3 (Pal 3 younger) 
if their abundances are the same." Fig. 12 clearly shows 
that such a larger TO age difference is not favored by the 
present models under the same assumption as adopted by 
Stetson ct al., namely: identical abundances. Age dif- 



ferences for this pair in the 2 Gyr range might still be 
reconciled with the HB morphologies of the studied clus- 
ters, provided other "second parameters" vary in such a 
way as to "compensate" the effects of such relatively large 
age differences upon the predicted HB morphologies. The 
a-element abundance is one such plausible candidate; in- 
deed, [a/Fe] has been observed to be lower in the "young," 
loose GCs Pal 12 and Ruprecht 106 (Brown, Wallerstein, 
& Zucker 1997) than in the vast majority of Galactic GCs 
(Carney 1996 and references therein). However, a lower 
[a/Fe] for Pal 3 would drive a bluer HB type for the clus- 
ter, thus decreasing even further the (already small) age 
difference needed to account for the difference in HB mor- 
phology between M3 and Pal 3. Hence a lower [a/Fe] for 
Pal 3 would work in the opposite sense, from the stand- 
point of HB models, to what would be needed to produce 
agreement between the present results and a TO age dif- 



E 
_o 

O- 

I 



1.50 



1.25 - 



1.00 



0.75 - 



0.00 



1.75 



1.25 - 



1.00 - 



n r- 



0.00 



— ■ ■ i — 



Eq. (A1) 

Eq. (A3) 
Eq. (A4) 
Eq. (A5J 




a) Reference: redder HB of M3 (r > 21 0") 
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b) Reference: bluer HB of M3 (r < 50") 
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Fig. 12. — The difference in age AAge(M3 — Pal 3) (in Gyr), derived for the several indicated mass loss formulae (see the Appendix in Catc- 
lan 2000), is plotted as a function of the M3 age (also in Gyr). In the lower panel, the best-fitting M3 simulations for the innermost (r < 50") 
radial range are considered (bluest HB type), whereas in the upper panel the best-fitting M3 simulations for the outermost (r < 210") radial 
range arc employed (reddest HB type). 



ference of about 2 Gyr. 

Of course, it may also be that future observations will 
reveal that the [Fe/H] values of Pal 3 and M3 differ, thus 
removing these clusters from the "bona fide second pa- 
rameter pair" category which they currently seem to oc- 
cupy. Empirically, one can see, from an HB morphology- 



mctallicity diagram such as the one presented in Fig. 1 
of Catelan & de Freitas Pachcco (1993), that an increase 
in [Fc/H] by 0.2 dex (i.e., within the Armandroff et al. 
1992 derived uncertainty) could move the Pal 3 HB to the 
red, compared to that of M3, by approximately the right 
amount, even if no age difference is present. (We have, of 
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course, ignored, for the sake of argument, both the inter- 
nal second parameter in M3 and the intrinsically narrower 
mass distribution of Pal 3.) Conversely, it is conceivable 
that, if Pal 3 turns out to be more metal-poor than M3 
by 0.2 dex, a larger TO age difference with respect to this 
cluster will be required in order to reproduce its red HB. 
Unfortunately, our imperfect (at best) knowledge of the 
variation of RGB mass loss with metallicity prevents us 
from carrying out a more meaningful analysis of the ef- 
fects of metallicity variations upon the HB morphology of 
Pal 3; again, the importance of improving our knowledge 
of mass loss in first-ascent giants cannot be overempha- 
sized. 

5. CONCLUSIONS AND DISCUSSION 

In the present paper, we have carried out a detailed anal- 
ysis of the "second-parameter pair" of globular clusters M3 
and Palomar 3, based both on observational data and the- 
oretical models. We have shown that the HB morphology 
of M3 is significantly bluer in its inner regions (observed 
with HST) than in the cluster outskirts (observed from 
the ground), and have presented Monte Carlo simulations 
which support the view that this difference in HB type is 
unlikely to be due to statistical fluctuations. This result 
leads to the conclusion that some "local" second parame- 
ter^) is acting between the inner and outer parts of M3. 

Our Monte Carlo simulations also suggest that the dis- 
persion in mass of the Pal 3 HB is very small (consistent 
with zero) , which is unlikely to be due to a statistical fluc- 
tuation in spite of the small sample size in the cluster 
CMD. In this regard, the difference in HB mass dispersion 
between M3 and Pal 3 is clearly found to be statistically 
significant. 

Ignoring the internal second parameter in M3 and the 
narrow mass distribution in Pal 3, we can do a traditional 
analysis searching for a global second parameter acting be- 
tween the two clusters. For example, we can search for an 
age difference based on the mean color or HB mass of each 
cluster. In this case, under the assumption that the two 
clusters have the same metallicity, the relative HB types of 
M3 and Pal 3 can easily be accounted for by a fairly small 
difference in age between them, of order 0.5-1 Gyr — which 
is in good agreement with the relative age measurement, 
based on the clusters' turnoffs, recently presented by Van- 
dcnBcrg (2000). Hence age could be the primary "global" 
second parameter for the pair M3-Pal 3. 

Even though there may be a small age difference be- 
tween M3 and Pal 3 (Stetson et al. 1999; VandenBerg 
2000), the dramatically different HB morphologies obvi- 
ous from a cursory examination of the CMDs appears to 
arise from another factor — the small mass dispersion of 
the Pal 3 HB. For some reason the star-to-star differences 
that lead to gm ~ 0.02 — 0.03 M Q in most GCs do not 
seem to occur in Pal 3. What possible mass loss-driving 



stellar parameter could have such a behavior? 

M3's internal second parameter adds additional com- 
plexity to the problem of identifying a global second pa- 
rameter acting between M3 and any other cluster. It is 
unlikely that age is the internal second parameter because 
over the > 10 Gyr age of the cluster dynamical mixing 
should have eliminated any initial radial separation of two 
age groups. On the other hand dynamical separation of 
stars which have undergone differential mass loss is un- 
likely to occur on a time scale comparable to or less than 
the HB lifetime. Even if it did, the settling of the more 
massive red HB stars to the center would produce a result 
opposite that observed. The most plausible internal sec- 
ond parameter is that some environmental factor leads to 
enhanced mass loss for the inner stars. One can at least 
imagine that stellar encounters could affect rotation rates 
which could in turn affect mass loss. If this is the case other 
clusters with HBs similar to M3, e.g., M5, might show an 
internal second parameter; it would be interesting to inves- 
tigate this possibility further using a combination of HST 
and ground-based photometry (as done in this paper for 
M3). Perhaps the density of the stellar medium, at least 
in this case and in Pal 3's, will turn out to be the culprit, 
provided it can be shown that RGB mass loss and the den- 
sity of stars in globular clusters are somehow related (with 
higher density leading to an increase in the overall amount 
of mass loss on the RGB and in the dispersion around the 
mean value). At least in the case of M3, another intrigu- 
ing possibility is that the detected radial variation in HB 
morphology is somehow related to the peculiar radial vari- 
ations observed in the blue stragglers populations in the 
same cluster (Ferraro et al. 1997b). Again we find our- 
selves in the position where a detailed analysis of a second 
parameter pair does not lead to a simple answer. 
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